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Abstract 

The cosmology in the Hubble expansion era of the Horava-Witten 
M-theory compactified on a Calabi-Yau threefold is studied in the 
reduction to five-dimensions where the effects of the Calabi-Yau man- 
ifold are summarized by the volume modulus, and all perturbative 
potentials are included. Matter on the branes are treated as first or- 
der perturbations of the static vacuum solution, and all equations in 
the bulk and all boundary conditions on both end branes are imposed. 
It is found that for a static volume modulus and a static fifth dimen- 
sion, y, one can recover the four dimensional Robertson Friedmann 
Walker cosmology for relativistic matter on the branes, but not for 
non-relativistic matter. In this case, the Hubble parameter H becomes 
independent of y to first order in matter density. This result holds also 
when an arbitrary number of 5-branes are included in the bulk. The 
five dimensional Horava-Witten model is compared with the Randall 
Sundrum phenomenology with a scalar field in the bulk where a bulk 
and brane potential are used so that the vacuum solutions can be 
rigorously obtained. (In the Appendix, the difficulty of obtaining ap- 
proximate vacuum solutions for other potentials is discussed.) In this 
case non-relativistic matter is accommodated by allowing the distance 
between the branes to vary. It is suggested that non-perturbative po- 
tentials for the vacuum solution of Horava-Witten theory are needed 
to remove the inconsistency that non-relativistic matter creates. 



1 Introduction 



Over the past several years there has been much interest in examining cosmol- 
ogy in five dimensions to see if the conventional four- dimensional Robertson 
Walker Friedman(RWF) cosmology can be recovered. Theoretical motivation 
for such attempts can be found within the framework of the Horava-Witten 
M-theory the strong coupling limit of the heterotic string theory. 

In Horava-Witten (HW) theory, space is eleven-dimensional (IID) bounded 
by two ten-dimensional orbifold planes with a Zi refiection symmetry in 
the eleventh dimension. In the lowest approximation Mio = M4xXx5'^/Z2 
where X is a compact Calabi-Yau (C-Y) threefold. Eleven dimensional super- 
gravity resides in the bulk. The construction of a consistent theory depends 
on a set of interlocking constraints due to anomaly cancelation, gauge in- 
variance, and local supersymmetry invariance. Thus there must be E'^ gauge 
interactions with chiral multiplets on each orbifold plane (one being the phys- 
ical world, the other the hidden sector) and a consistent theory exists only if 
their exists a relation between the ten-dimensional gauge coupling constant 
and the eleven-dimensional Planck mass. Assuming the compactification 
mass is the GUT mass, Mc this latter relation then shows that the eleventh 
dimension is (9(10) times larger than the C-Y scale. Dimensional reduction 
then leads naturally to a five dimensional world M^^S^ / Z2, with the residual 
effects of the C-Y manifold being described by their moduli. The details of 
this formulation have been extensively developed in [5]. 

Phenomenological analysis of five- dimensional cosmology was stimulated 
by the work of Binetruy, Deffayet, and Langlois [HI, and subsequently by the 
Randall-Sundrum (RS) model [71IH]. In the RSI model, where space is five 
dimensional with five- dimensional gravity in the bulk and physical matter on 
two four dimensional bounding planes, it has been argued that the standard 
four-dimensional RFW cosmology is obtained only if the fifth dimension is 
stabilized [0] (though counter arguments have been given in JO]). In order to 
achieve this naturally, it has been suggested that one phenomenologically add 
a scalar field in the bulk with appropriate bulk and brane potentials [TT| IT^ . 
Vacuum solutions appropriate for these potentials are then obtained and 
solutions are discussed with brane matter viewed as perturbations on the 
vacuum [131 Ej. 

In this paper we examine the five-dimensional reduction of HW theory 
in the simplest approximation of keeping only the volume moduli of the 
C-Y space, and consider under what circumstances one can reproduce the 
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standard RWF cosmology. When no matter is present on the branes, a 
vacuum solution of the gravitational equations has already been given in [H] . 
This vacuum solution automatically fixes the bulk and brane cosmological 
constants (without any fine tuning) so that when matter is added to the 
branes the net cosmological constant is correctly zero. Treating the matter 
as a pertubation on the vacuum, we find that the standard cosmology can 
be recovered for relativistic matter for a static solution but not for non- 
relativistic matter. We further show that this result remains even when an 
arbitrary number of five-branes are added in the bulk. We then examine 
the RS phenomenology, which closely resembles the 5D reduction of HW M- 
theory. Vacuum solutions for this system are difficult to find for an arbitrarily 
chosen bulk potential. Further, the difficulty arising from non-relativistic 
matter in the HW theory is evident only when the boundary conditions on 
both orbifold planes are rigorously imposed. Therefore, we examine the class 
of models presented in which can produce exact vacuum solutions for a 
specific type of potential. We show that the inconsistencies that arise in the 
HW theory due to the presence of non-relativistic matter can be avoided due 
to the introduction of ad hoc potentials not found in the HW M-theory with 
a subsequent fine tuning needed to set the brane cosmological constant to 
zero. (In contrast, HW theory sets the cosmological constant to zero without 
any fine tuning.) The solution is then static only to lowest order, the brane 
separations varying in time due to the presence of non-relativistic matter. In 
an appendix we examine the RS model with the bulk potential given in ^3] , 
and discuss the difficulties in finding the vacuum metric for this case. 

In Section 2 we briefly review the Horava-Witten model and its 5D re- 
duction, and write down the field equations and boundary conditions. In 
Section 3 these equations are solved in the presence of matter, and we show 
that non-relativistic matter is excluded unless one fine tunes the matter on 
the two branes. In Section 4 we introduce 5-branes in the bulk (the only ad- 
ditional freedom allowed in HW theory) which reduce to 3-branes in the 5D 
reduction, and show the same problem remains (unless the non-relativistic 
matter on the end branes and on the bulk 5-branes is fine tuned). In Section 
5 we relax the static condition imposed in earlier sections and find that the 
fine-tuning of the matter on the branes might be alleviated due to the non- 
trivial bulk dependence of the time derivative of the scalar field. In Section 6 
we discuss the RSI model and show how it differs from HW theory. Conclu- 
sions are given in Section 7. An appendix discusses the choice of potentials in 
RSI used in ^|. It is shown that the approximate solution for the vacuum 
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metric used there does not give rise to a hierarchy when all the boundary 
conditions on both branes are imposed. 

Previous work on HW cosmology has been given in P^IT^ ITTj. However, 
the first two papers do not impose all the boundary conditions, and hence 
do not see the difficulties found here. The last paper is concerned with the 
inffationary era rather that the Hubble expansion era being discussed here. 
A very general analysis for an arbitrary model was given in ^H], but the 
authors do not seem to have noticed the difficulty discussed here. Within 
the M-theory framework, there have been several papers suggesting that 
stabilization of moduli can be achieved by turning on fiuxes fSll2nil^- The 
first two give rise to large negative cosmological constants, while the last 
to a large positive cosmological constant. How to reduce these constants to 
their physics values in a natural way remains (as well as how to analyse the 
presence of matter). 

2 Horava-Witten M-theory 

In this section we review the basic formulae of Horava-Witten M-theory. Ho- 
rava and Witten ^1121 showed that the low energy limit of eleven dimensional 
supergravity on the orbifold R^^^S^ / Z2 is dual to the strong coupling limit 
of the ten dimensional x heterotic string. This theory is reduced to five 
dimensions via compactification on a Calabi-Yau threefold of volume V. In 
order to cancel the gauge and gravitational anomalies that arise, an E^ gauge 
group is required to reside on each of the two lOD planes at the orbifold fixed 
points (x^^=0 and vrp), and a relation between the Yang-Mills gauge coupling 
constant and the eleven- dimensional gravitational coupling constant is pro- 
duced . This is remarkable in that it explains the discrepancy between the 
GUT scale Mg— 3x lO^^GeV and the four- dimensional Planck mass, since it 
is the eleven- dimensional Planck mass, which is ~ M^, that is the funda- 
mental scale. Upon reduction to four dimensions and using the values of 
the four-dimensional Newton constant Gat, as well as acuT and M^, one 
finds that the orbifold radius is about ten times the compactification scale, 
rendering the theory effectively ffve-dimensional. 

Discarding the shape moduli and keeping only the volume modulus, V , 
of the C-Y threefold, the bosonic part of the reduced five- dimensional La- 
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grangian takes the following form jH] 
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where K5 is the five-dimensional Newton constant, R is the five- dimensional 
Ricci scalar, ac is the GUT scale coupling, Fj^^ are the gauge field strengths 
on the boundary orbifolds, C" are complex scalars of chiral matter, W is the 
superpotential, and the last term represents the D term of the gauge theory 
on the branes. The parameter a, which is (9(10^^GeV), fixes the bulk and 
brane cosmological constants. It is defined as [22] 
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where R is the Calabi-Yau curvature constructed from the metric gai, and 
is the Kahler form 

^ab = i9ah (2) 

with holomorphic and anti-holomorphic indices a and h. 
From this action and the ansatz 



ds^ = a{t, yfdx^dx^ - n{t, yfdt^ + yfd'if, 



(3) 



where y=x^^ is the coordinate of the fifth-dimension which extends from 
y=0 to y=vrp, and t is time, one can obtain the five- dimensional Einstein 
field equations 
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where i = 1, 2 corresponds to the fixed points at y=0, 7rp, prime and dot 
denote derivatives with respect to y and t respectively, and V = e*^ where 
<p is the breathing modulus of the Calabi-Yau. The non-relativistic matter 
density on the i'th orbifold is pinr-, the relativistic matter is pir, and is the 
pressure. 

The 5-functions in Eqs.(4) and (5) imply boundary conditions at the 
orbifolds y=0, vrp given by 
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where M5 is the five-dimensional Planck mass and pi and pi are the total 
matter density and pressure on the two orbifolds. Thus the bulk cosmological 
constant a, and the brane cosmological constant ctj are naturally correlated 
without any fine tuning. 

In addition to the Einstein field equations, one can derive field equations 
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and boundary conditions for the breathing modulus from the action: 




0. 







(11) 



(12) 



This differs from the resuh of Ref.fTH] in that non-relativistic matter has 
been included in the boundary condition for cj) as it should since is coupled 
to the gauge fields on the branes at y=0 and vrp by the factor V^~^=e~'^ in 
the superpotential term in Eq.(l). (Note that also couples to gauge fields 
in Eq.(l) with the factor V=e'^, but the coefficient F^^F^^ vanishes for the 
radiation fields). 

3 Solution of the 5D equations 

We now proceed to solve the field equations using a perturbative expansion 
in powers of matter on the branes. This expansion is allowed since in the 
Hubble era the matter density is very small compared to the cosmological 
constants in the bulk and on the brane which are of GUT size. We start with 
the vacuum solutions and then include matter on the branes as a higher order 
correction. The vacuum solution which fully solves the bulk and boundary 
equations, preserves Poincaire invariance, and breaks 4 of the 8 supersymme- 
tries (appropriate for getting N=l supergravity when one descends to four 
dimensions) was given in [5] as 



Here ho is a constant that is arbitrary due to the flat directions of the potential 
and f is given by 



where c is a constant. To first order in p, the vacuum solutions are perturbed 
to take the following form: 



a{y) = ; n{y) = ; b{y) = Kf ; V{y) = Kf. 



(13) 



f{y) = c + a\y 
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n{y) = + (16) 
h{y) = h^fil + 5h{y,t)) (17) 
V{y) = h^f{l + 5V{y,t)) (18) 



It is convenient to introduce tlie notation 



Aa' = 5a' + ^(51^-^56 (19) 
An' = 5n' + — - —56 (20) 

IW' = 5V' + y5V-^5h (21) 



along with the definition of the Hubble constant 



a 



H = -. (22) 
a 

The significance of the combinations of Eqs.( 19-21) is that they are invariant 

under a first order coordinate change in the y coordinate: y = J + 

To first order in p, the Einstein equations Gu, G^k, Gyy, and the field 
equation for the breathing modulus become the following: 

Aa" = blf (h' + hI- ^ hlfA, (23) 



An" + 2Aa" = hlf [?,H^ + 2H + 2H^- + ^- + ^ blfA^ 



(24) 



3Aa' + An' - Ay' = ^ 4//2 + 2// - ^ = -^A^, (25) 

ay 6 y a 

Ay" + y (An' + 3Aa' - AV) = hlf + 3i/0 + = hlfA^ (26) 

As we will show in Section 5, the Gty equation is of higher order and will be 

discussed later concerning the possibility of inchiding y-dcpcndcncc in the 
Hubble constant. Inserting the metric ansatz into the boundary equations 
(8), (9), and (12) yields 
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A«i = (-1)'" (2ft, + 3p.,) + ^iiJ:P.nr^ (28) 

where fi are the values of f(y) on the branes: 

/i = c ; f2 = c + a2np (30) 
It will also be helpful to use the combination of Eq.(27) and Eq.(28) that 

isolates Pinr- 

(3Aa' + An'),, = ^"^^'^^t^'"-^ ^^^^ 

where we have used Pir—Pir/^- The significance of the above results is that 
both the field equations and the boundary conditions can be expressed in 
terms of the y- invariant combinations of Eqs.( 19-21). 

Wc arc now ready to examine the solution of the bulk equations, impose 
the boundary conditions on them, and check their consistency First of all 
we notice that the equations Gu and Gkk are easily solved by integration 
with respect to y. However, it will be more convenient to change integration 
variables from y to f and to use the combination Gu + Gkk instead of Gkk- 
One finds 

Aa' = I ( /^^ df'f'A, + ci j (32) 
3Aa' + An' = I ^J^^ dff {A^ + A^) + ci + (33) 
and imposing the boundary conditions Eqs.(27) and (28) one obtains 

c, = -^((.;/,V„ + iA.„) ; A^il^ (34) 

f^' dfTA, = ^ (bl (p,,/? + + + (35) 

Cl + C2 = --^Plnr (36) 
0/1 

dff {A, +A,) = ^ f Ipi,, + lp2nrl • (37) 
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Eqs.(35) and (37), arising from the boundary conditions at the distant 
brane y=7rp, thus produce constraints on the time derivatives (which enter in 
the Ai defined in Eqs. (23-26)) of the metric and 0. We next integrate the field 
equation for the breathing modulus, Eq.(26), using Eq.(25) to eliminate 3Aa' 
+ An' -AV. Combined with the boundary condition Eq.(29) one obtains 

AV'=^^^J^yf'f"{A,-3As) + cs^ ; C3 = ^pin. (38) 



ot \Jfi J h 

j^^'^ dff (^4 - 3A3) = -A l^j^p^nr + ^P2nr] • (39) 

g equation to be satisfied is 
(38) back into Eq.(25) gives the constraint 

'y^^ dff (^1 + ^2 + 3^3 - ^4) + Ci + C2 - C3 ) = ^^3 (40) 



Jl J2 

The remaining equation to be satisfied is Eq.(25). Inserting Eqs. (33) and 



This constraint is a strong one as it must hold for all y. 

We now examine the consistency of this system. We consider here the 
static case where (j) — — b. Here ^4=0 and = 741+742. Thus multiplying 
Eq.(37) by 3 and adding to Eq.(39) gives 

= (^Y^P^nr + ^P2nr^ (41) 

Thus a consistent solution without fine tuning requires (when (p = — b) 

Pnr = (42) 

i.e. only relativistic matter is consistent with Horava-Witten cosmological 
equations. However, in addition to Eq.(42) one must also make sure that 
the constraint Eq.(40) is satisfied. In Section 5 we will show that the Gty 
equation implies H' is 0{p^^^) for the static case, and hence to 0{p) that we 
are calculating one can consider and H to be independent of y. Hence we 
note that for the static case Eqs. (35) and (37) correctly reduce to the RFW 
cosmology equations for relativistic matter with Gn defined in terms of A 
and fi-. 

= y GiV {P[r + P'2r) ; = ^ {jf^Tt) ^^^^ 
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H'^ + 2H = (44) 

(where the rescaled = hoffpir are the mass densities as seen locally on 
the orbifold 3-branes). Eqs.(43) and (44) just incorporate the 4D relativistic 
matter equation of continuity: p = -3H{p + p) = -4 Hp. Since Ai + A2 
= A3, Eq.(40) is then identically satisfied as a consequence of Eqs.(42) and 
(44). 

In summary we note that it is the boundary conditions on the distant 
brane at y = irp, Eqs.(37) and (39), that produces the constraint Eq.(42) 
on non-relativistic matter (which is why earlier analyses have not seen this). 
However, a satisfactory FRW cosmology does result for relativistic matter, 
with the brane cosmological constant naturally vanishing with no fine tuning 
required. 



4 Inclusion of 5-Branes 

We have shown that in the static case (0, 6 = 0) the system of bulk equations 
with their boundary conditions imposed is inconsistent when non-relativistic 
matter is included in the system. Therefore we would like to examine other 
situations that might lead to a consistent solution when all types of matter are 
present. The only additional generalization available in the Horava-Witten 
theory is to include a set of 5-branes in the bulk transverse to the orbifold 
direction 01231 • We follow here the analysis of of a single 5-brane residing 
at an arbitrary position y=Y in the bulk (which can easily be generalized to 
an arbitrary number of 5-branes). The fields that live on the 5-brane include 
an N=l chiral multiplet and N=l gauge multiplets but no superpotential. 
One must generalize the function f(y) and the definition of a to be (for 
0<y<7rp) 

f{y) = c + h{y) ; h{y) = -a,y + {~a^y + a,Y) e{y-Y) (45) 

h\y) = a{y) = -«i - a.O {y - Y) (46) 
a{y = Y)^as = ^ (47) 

with the cohomology condition 

ai + a2 + a5 = (48) 



10 



but otherwise the vacuum solution has the same form as Eq.(13). Note that 
f(y) is continuous whereas Q;(y) = ^'(y) is not 



\ -CKi - CKs = q;2 y < y < 7rp ^ ' 

Wc use the same ansatz as before for the metric: Eqs. (14), (15), (16), (17), 
and (18). However, the Einstein equations are altered due to the fact that 
a is no longer a constant; thus a' no longer vanishes. We also must make 
some assumptions about the matter content of the five-brane in the bulk. 
We know that in order to give rise to the Big Bang at the end of inflation 
that the inflaton couples to matter on the physical orbifold at y=0 and we 
have assumed that it also couples to matter on the orbifold at y=7rp. There 
is no a priori reason to believe that it also couples to any matter fields on 
the 5-brane in the bulk. However, we will make the assumption that it does 
couple to the five-brane matter with the same strength V as for the two 
orbifolds (this does not effect the general conclusions of this section). 

Now we would like to solve the Einstein equations in the presence of this 
5-brane. The calculation is very similar to that done in the previous section 
with the modifications that i runs from 1 to 3 and a' terms must now be 
included. For example in the Gu equation we still have the definition 

Aa' = 6a' + —6V - —M. (50) 
2/2/ ^ > 

However Aa" is 

Aa" = 5a" + ^ - ^ - ^ i^b' - SV) ' ^ (Sb - SV) (51) 
and after making use of the ansatz, the Gtt equation becomes 
Aa" = (p,rblf^ + -ip3n.) S{y-Y) + hlfA^ 



(52) 



where the 5-function term arises from a y-derivative of a, and Ai is given in 
Eq.(23). The key point in this relation is the unexpected result that the 
terms cancel, not only at the vacuum order (which was already verified in 
showing Eq.(13) with f(y) of Eq.(45) is correctly the vacuum solution), but 
also at higher orders. This also applies to the other field equations. 
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Since there is a discontinuity in a when crossing the five-brane, we solve 

the Gtt equation in two domains: 0<y<Y and Y<y<7rp and then match them 
across y=Y. We define the gauge covariant combinations in these regions as 

Aa[ = 5a' + ^ {Sb - 5V) ; 0<y<Y (53) 

Aa2 = 5a - {5b - SV) ; Y <y< np. (54) 

These quantities obey the boundary conditions 

Aa; (y = 0) = - ^ (^p,rblf! + jp,^r^ (55) 
Aa'^iy = ^P) = ^ (^P2rblf2 + ^Psnrj • (56) 
We now solve the Gu equation in these two domains to obtain 

Aa[ = £ dy'blfA, - ^ [p.rblfl + j^p,n)j (57) 

A4 = - fj dy'blfA, + ^ (^p,rblfl + j^p2n}j . (58) 



The discontinuity across the five-brane required by Eq.(52) gives 
j^^^y'Aa" = (A4 - Aa\)^^y = --^ (^P^rhlfl + j^p^nr^ 
Thus, subtracting (57) from (58) leads to 



dy'blfA, = E ^ (pM! + -^AnrJ (60) 

If we make the assumption as before that = = 6 this equation will give 
the RFW relation for the Hubble parameter but now with matter from three 
separate branes included. However,the situation with no matter included on 
the 5-brane in the bulk does not reduce to the case with only two branes 
since a and therefore the function f(y) are modified from the case where only 
two branes were present. Therefore, the Hubble law (more specifically the 
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Newton constant, Gjv) is affected by the presence of the additional five brane 
even if it is empty. 

The other field equations can be solved in a manner similar to that de- 
scribed above for the Gu equation. Namely, we first modify the equations 
with the inclusion of terms involving the y-derivative of a, look at the equa- 
tions separately in the regions, 0< y < Y and Y < y < 7rp, and then match 
them across the five brane. The result for the (j) equation in the bulk is 

Ay" + J (An' + 3Aa' - AV^') = hlfA, + -^^P^nrS {y - Y) (61) 



which leads to 



■Kp ^' 1 

dy'hlf {A, - 3A3) = - E J^P^nr. (62) 



The Gkk equation gives 

j dy'hlf {A, + A,) ^Y.-^^ Pinr. (63) 

The Gyy equation remains unchanged from Eq.(25) since there are no matter 
sources present but where now the gauge invariant combinations Aa, An, 
and /SV include the new definitions of a{y) and f(y). We can now proceed 
to check the consistency using the same steps that led to Eq.(41). The new 
relation is 

-rPlnr + -rP^nr + -y P^nr = (64) 
\/l h h / 

Once again we sec that introducing non-relativistic matter into the system 
results in an inconsistency if there is no fine tuning of the matter on the 
different branes. 



5 y-Dependence and Non-Relativistic Mat- 
ter 

Assuming the static case for and b we have shown that the system does 
not admit consistent solutions in the presence of non-relativistic matter. We 
would now like to relax this assumption and then examine the system. We 
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consider here the case of no 5-branes present. To see what constraints are 
put on our assumptions, let us first look at two separate ways of evaluating 
the y-derivative of H and compare this with the Gty equation. The definition 
of H is 



H = -. (65) 
a 



Therefore 



d' a' a 



H'^ (66) 

a a a 



and using the explicit form of a{y, t) in Eq.(15) we find 

b! ( a 



Alternately 
which gives 



o! a 



- = - + 5a' (68) 
a 2/ 



^= f^ + 5a'ji/ + 5a'. (69) 



Comparing Eq.(67) with Eq.(69) and using the fact that in the static case 
Aa' = (5a' we find 

Aa' = E'. (70) 
In the static case the Gty equation is given by 

Ad' = (An' - Aa') E. (71) 

The right hand side is seen to be of order p^l'^ which shows from Eq.(70) that 
E' is also of order p^/^. Thus the static case requires H to be independent of 
y to first order in p. 

Let us next examine the situation when we let and b depend on time. 
The Einstein equations Gtt, Gkk + Gu-, Gyy, and Gyy plus the 4> equation of 
motion (Eqs.(35),(37),(39),(40)) give 

j;;m - ^,±}iP.f!^f) (72) 

/^df'nAr + A,) = (73) 
14 



2 

Ji ,-1 Ji 



^'rf/'r(3A3-A4) = (74) 



^ dfnA, + A, + 3A, - A,) = f'A,{y) + ^-^^ (75) 

fl Ji 



In particular evaluating Eq.(75) at y = ?/i gives 

Mvi) = (76) 

However, considering the combination Eq.(73) + Eq.(74) - Eq.(75), where 
Eq.(75) is evaluated at y = y2, yields 



Mv.) = ^ (77) 



and therefore we see that A3 has non-trivial y-dependence. Since ^3 only 
depends on H and this implies that depends on y. With the static 
condition relaxed we can now see that the Einstein equations contain integrals 
over functions whose y-dependence is not known. Therefore relations such 
as Eq.(41) are no longer valid. Eq.(41) now becomes 

J^^' df (^9H^- + 3^ + 0^ + + ^0 + 3H<pj = (J-p,^r + j-p2n?j (78) 

and we see that non-relativistic matter is now related to an integral over an 
unknown function of y. Without knowing the exact form of the integrand 
it is difficult to determine if there is any constraint on the non-relativistic 
matter with the static constraint relaxed. 



6 Randall- Sundrum Model 

The Randall Sundrum Model (RSI) |7j is a phenomenological five dimen- 
sional model with the fifth dimension bounded by two 3-branes, and thus 
it resembles the 5D reduction of Horava-Witten M-theory. It is of interest 
then to compare the two to see what differences exist. A problem arose in 
the earlier versions of RSI when it was realised that the correct Hubble era 
cosmology could not be reproduced, and it was then suggested that one could 
obtain the RWF cosmology if one were to stabilize the fifth dimension [D]. 
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A technique for this was proposed by Goldberger and Wise [TT] by adding a 
scalar field in the bulk with appropriate bulk and brane potentials, and this 
idea was further analyzed by Cline and Firouzjahi [T^. However, neither 
Refs. or ^3] included the dynamics of the gravitational field in the bulk. 

In general it is very difficult to find vacuum solutions for the Einstein- 
scalar field equations with arbitrary brane and bulk potentials. Further, we 
have seen in Sec. 3 how important it is to have a rigorous vacuum solution 
since it is the careful imposition of boundary conditions on both branes 
that produces tension in the system. However, a special class of solutions 
for the vacuum metric obeying all boundary conditions on the branes was 
constructed by de Wolfe et. al J2] where the brane and bulk potentials are 
related to a single function of the scalar field. Subsequently, in Ref . [TB] . 
matter on the branes was added and treated perturbatively with respect to 
a vacuum solution. However, an explicit vacuum solution for their choice 
of potentials was not obtained, and the difficulty in obtaining a solution is 
discussed in the Appendix. Here we will discuss the case of [12] with matter 
treated as a perturbative addition on the branes , since this case treats the 
vacuum solution rigorously. 

In order to make contact with the notation found in [T2j we let V^(0) be 
the bulk potential and Aj(0), i=l,2 be the potentials on the two branes at 
j=yi and J=y2 respectively. The metric is given by 



(79) 



with the perturbative expansions 



N{t,y) = Ao{y) + 6N{t,y) 

A{t,y) = Ao{y) + 6A{t,y) 

b{t,y) = l + 5b{t,y) 

(t){t,y) = <Po{y) + H{t,y). 



(80) 
(81) 
(82) 
(83) 




(84) 
(85) 



(86) 
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where primes denote ^ except onV{(j)) where it represents The boundary 
conditions are given by 



A' 



y=yi 



y=yi 



' 2 d(h. 



y=yi 



y=yi 



(87) 
(88) 



where i = 1,2 corresponds to the brane locations at y = yi, y2- 
It was shown in [12] that if V{(j)o) has the form 



1 fdW{<j)o) 



(89) 



for some W{(j)o), then a solution of the bulk vacuum equations will also satisfy 



_ l dW{<Po) 
'^^ 2 90„ 

as long as the boundary conditions 



(90) 



y=yi 



y=yi 



y=yt 



y=yi 



(91) 

are also satisfied. Once Vr(0o), Ai, and A2 are chosen one is then left with a 
set of first order differential equations that can easily be solved. 

In the example constructed in [T2] W^(0o), Ai, and A2 were chosen to be 



m0o) = 

Ai(0o) = w{ct>o{yi)) + w\Uyimo-Uyi)) 

+7i(0o - 'Po{yi)f 
A2(0o) = -w{Uy2))-w'{Uy2mo-Uy2)) 

+72(0o - (po{y2)f 



which gives 



0o(l/) 

Ao{y) 



he 



-liy 



ao - f - \<Ple-'^' 



(92) 

(93) 
(94) 

(95) 
(96) 
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where Co, L, 71, 72, 0i, and 02 are all arbitrary parameters. The relationship 
between Aj(0o) and W{(f)o) fine tunes the net cosmological constant on the 
branes to zero. While there is no a priori motivation for the choice (other 
than the fact it fine tunes the cosmological constant to zero), it leads to 
simple analytic forms, Eqs.(95),(96), which can be treated easily. 
Next we define the quantities 



AN' 
AV 



5A' + 
SN' + 



-A'Jb 



L'2 



5b. 



(97) 

(98) 
(99) 



These quantities are invariant under a first order coordinate change in the y 
coordinate: y = j + S(j). 

To first order in the matter the Einstein equations Gu, Gkk, Gyy-i and the 
field equation for in the bulk become 

AA" + AA'^AA' = e-2^°(f)ij2 = (iqO) 
2AA" + 8A'^AA' + AN" + 4A'^AN' - e-^^^^^^ {SH^ + 2H) = A2 (101) 

A'„ {3AA' + AN') - ^AV = e-^^"^^) (2H^ + h) = A3 (102) 
AV" + AA'^AV + {3AA' + AN') = 0. (103) 
The boundary conditions are 

AA' = {-ly ^ (104) 
y=yi 3 

(3AA' + AN') = (-1)' ^ (105) 



y=yi 



AV 



y=yi 



= 50, ((-1)^+S^0; - 



(106) 



y=yi 



where 7, are the free parameters in the model of Eqs.(93),(94) i.e. 



1 92 A,, 



7i = o 



2 



(107) 



Immediately we sec a difference between HW theory and the RS phe- 
nomenology manifesting itself in the boundary conditions for AV. In HW 
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/■2/2 

/ 



only the invariant quantity AV appears in the boundary conditions whereas 
an additional 6(f) term appears in the RS model. One can easily check that 
in HW there is no free parameter analogous to 7^, and the M-theory choice 
of potential makes the term in parenthesis on the right-hand side of Eq.(106) 
zero. This distinction will allow one to avoid the constraint on matter found 
in the HW theory. 

We will now solve the system in the presence of arbitrary matter. The 
Gtt equation can be integrated to obtain 

AA' = e-^^°(^) r dy'e^^°^y'^H^ + e'^^-^^'^ci (108) 
Jyi 

and using the boundary conditions we get 

ci = -Ple^Myi) (109) 
3 

^yf^2A4y')H^ = Ple^Myi) + Ple^My2) (ng) 
'yi 3 3 

If H is independent of y one can see that the usual Friedman relation is 
recovered. Eq.(108) now becomes 

AA' = e-^^°(^) ( r dy't^^'^'^y'^H^ - ^e^^°(^^)) (111) 

In a similar manner we can solve the combination Gu + Gu to obtain 

dy'e^^°^y'\AH^ + 2H) = ^e^^°(^^) + ^e^-^-^^^) (112) 

3AA' + AN' = e-^^°(^) ( r dy'^^''^y'\AH^ + 2H) - ^e^^°(^^)) . (113) 

\Jy\ 3 / 

The remaining equation is Gyy which can be solved for AV or equivalently 
for Scpiy) in terms of 5h{ii) 

5m = 0; /' dy'6b{y') + 0^ /' dy'^ {3AA' + AN') 
Jyi Jyi (Po 

-^'^£dy'^^^ + h)+ 5^)- (114) 
The boundary conditions on (f){y) are^ 



^We note the limit 7^ ^ 00 is the stiff potential limit of 
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One can then substitute the known expressions for A'^, 0^ (for the specific 
model of [12] these are given by Eqs.(95),(96)), and 2H'^ + H from Eq.(112) 
into Eq.(115). This then determines S(f){yi) in terms of pinr- However, there 
remains one additional relation involving 6(f){yi) for we may let y = 1/2 in 
Eq.(114). Eliminating S(j){yi) and 50(|/2) using Eq.(115), then Eq.(114) at y 
= y2 becomes a constraint on Sb{y): 



~L dy'^^i'i^A' + ^N'). (116) 



Note that the right hand side of Eq.(116) depends only on the vacuum metric 
and Pinr which can be seen upon substitution of Eqs. (112), (113), (115) 

'"dy'5h{y') = 
yi 

K{y2) „ _ p-2Ao(j/2) p-l V Pina.p'^Aoiyi) 



+ 



2 



P1J9 /o~2ylo{y') ^ , 



yi 

Jyi 110/- i\ 



<P'^{y') 3 



where 



dy'e'^''^y'\ (118) 

'j/i 

Therefore a consistent solution for arbitrary non-relativistic matter is ob- 
tained. (It should be noted that although these equations hold in general, 
explicit values for S(f){yi) can only be calculated if A'^ and 0^ can be deter- 
mined, which will be dependent on the choice of the bulk potential.) 

We see now the meaning of the result that the 6(j) boundary condition de- 
pends on both the coordinate invariant combination AV and on 6(f)i (rather 
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than just on AV as in HW). Instead of piitting a constraint on pj„r, this de- 
termines the integral of 6b{y) in Eq.(116) which is just the change of distance 
between the branes due to the presence of non-relativistic matter. This is 
possible with the phenomenological potentials of the RS model, but not in 
HW where the theory determines the potentials to automatically cancel the 
cosmological constant. Since pinr decreases as l/a(i)^, the brane separation 
is actually time dependent at higher order in RS so the distance between 
the branes cannot be fixed, and so the RS model is also non-static. Here, 
however, the time variation of the invariant distance between the branes sets 
in at 0{p^ll) since p ~ Hp = 0{p^^^) while in HW the time variation sets 
in at 0{p). 

The above illustrates the difference between the phenomenology of RS 
and the theory of HW. In the RS model one is free to add arbitrary bulk 
and brane potentials (characterised here by 7^) for the scalar field while in 
HW the couplings of the volume modulus V are determined by the theory 
and one is not free to include ad hoc potentials. In our analysis of HW 
theory we have used all potentials that arise perturbatively. There are non- 
perturbative potentials in HW that have not been included that might relax 
the tension that non-relativistic matter produces. If this is the case, one 
would also have to modify the vacuum solution to take into account the 
additional interactions. 

7 Conclusions 

We have studied here the Hubble era cosmology in Horava-Witten theory. 
After compactification of the IID space on a Calabi-Yau threefold, the system 
reduces to a 5D theory, the fifth dimension, y, bounded by two 3-branes with 
gravity and a scalar field (representing the volume modulus) in the bulk, and 
gauge and chiral matter on the 3-branes. The field equations were solved in 
the bulk and the boundary conditions on both 3-branes were imposed. We 
have shown that for the static solution, the standard RWF cosmology arises 
for relativistic matter on the branes, but the field equations cannot allow 
non-relativistic matter. This result arises from the constraint of satisfying the 
boundary conditions on both 3-branes. The same result maintains if one adds 
5-branes in the bulk (the most general form of Horava-Witten M theory) . We 
have included all of the potentials that arise perturbatively in HW theory. 
However, there are non-perturbative potentials in the theory that have not 
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been included and may allow the introduction of non-perturbative matter. 
Once these potentials are included, the vacuum structure of the theory will 
be altered and one will need to find a new set of vacuum solutions. We also 
have not addressed the issue of moduli stabihzation that has been studied 
recently in the context of flux compactifications in [T^ I21j . 

We have also demonstrated the difference between the HW theory and the 
RS phenomenology. For the special class of potentials studied in [T2] the RS 
model can accommodate arbitrary matter on the branes and thereby repro- 
duce the RWF cosmology. This occurs due to the existence of free parameters 
in the brane potentials which can relax the constraint on non-relativistic mat- 
ter found for HW theory by allowing one to solve for the change in brane 
separation due to the presence of matter rather than putting a constraint on 
the matter itself. In HW theory one is not allowed to introduce such free 
parameters as the form of the brane potentials are fully determined by the 
consistency conditions of the theory. 

In the Appendix we have analysed the RS model discussed in [12] which 
includes a scalar potential that does not fall into the class of potentials stud- 
ied in ^21- The vacuum solutions in this case are given by a series expansion 
and we have shown that one cannot generate the desired solution to the hi- 
erarchy problem if one truncates the series after the first term as was done in 
(or if one includes the second term). The difficulty is that the vacuum 
solution must obey boundary conditions on both branes which precludes the 
hierarchy from developing. It is an interesting question of whether this diffi- 
culty is due to the truncation or whether formation of a hierarchy is sensitive 
to the choice of brane and bulk potentials. 
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A Appendix 

In section 6 we discussed the Randall- Sundrum model and showed that for a 
general class of bulk and brane potentials one can obtain vacuum solutions 
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and then find consistent solutions to the Einstein and scalar field equations for 
arbitrary matter introduced perturbatively on the branes. In this appendix 
we will analyse the situation for a scalar field obeying the potentials chosen 
in JHlji-G. in the bulk 

= (A.l) 

and on the boundaries Vi($) = mj($ — fj)^. Here rrii are the analogs of 7^ 
that arose in Sec 6. These potentials cannot be put into the form specified 
by Eq.(90) and therefore one will have to solve the second order differential 
equations to find the vacuum solutions. We will find that it is not easy to 
see how this choice of potentials leads to the desired hierarchy. 

Throughout this section we follow the notation of The metric is 

ds' = e-2^(*'^)rft2 - e-2^(*'^) ' Kt, y)W (A.2) 

i 

and the perturbative expansions are given by 

N{t,y) = Ao{y) + 6N{t,y) (A.3) 
A{t,y) = Ao{y)+5A{t,y) (A.4) 
b{t,y) = bo + Sb{t,y) (A.5) 

= ^oiy)+mt,y)- (A.6) 

The Einstein equations and the scalar field equation at vacuum order are 

\'2 _ U,I2 _27 2^2\ , .2,2 



= j^[K'-mX^i)+kX (A.7) 
a: = ^'h'^' (A.8) 

= 4A'^<!>'^ + mX^o (A.9) 

where is given by = 1/ where M is the 5D Planck scale and and k"^ 
are related to the bulk cosmological constant A by A = k^. Only two of 
these equations are independent since the third equation can be generated by 
taking the y-derivative of the first and then inserting the second. Therefore 
a solution of any two of these equations will necessarily satisfy the third. We 
look for solutions of the form 



Ao = ao + /3y+£a„e-2"'^^ (A.IO) 

n=l 

00 

$0 = E Cne-('"-^)"^ (A.ll) 



n=l 
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where a = ekbo- 

Ref.fTB] assumed that truncating the series at n=l represents a good ap- 
proximation. However, it is easy to see that when this truncation is inserted 
into Eqs.(A.7-A.9) one generates the higher terms of Eqs.(A.10,A.ll). We 
would hke to examine the question of whether the higher terms can be ignored 
thereby giving the resuhs found in ^B]- While the equations are non-linear, 
it is still possible to obtain recursion relations. Inserting Eqs.(A.lO) and 
(A. 11) into the vacuum equation (A. 8) generates relations between a„ and 
Cn, the first few of which are 



4ai = yc? (A. 12) 



16a2 = — 6ciC2 (A.13) 



36a3 = y (lOciCs + 9cy . (A.14) 

One can also see that f3 = kbo from Eq.(A.7). Inserting Eqs.(A.lO) and 
(A. 11) into Eq.(A.9) we find 



^ Cn («' (2n - 1)' + APa (2n - 1) - m^bl) e'^^"-!)?' (A. 15) 

oo 

= 8«2 J2 n(2m - l)a„c„e-(2"+'"-^)^ (A.16) 



n=l 



n,m=l 

For n=l this gives an equation for the coefficients of e~°^, 

,21,2 



a' + Aa/3 - mX = (A. 17) 



or 



6 = -2 + ^4 + ^ (A.18) 

where we have taken the positive root so that e > 0.(This is identical to the 
result found in .13 .) For n = 2, after using the result of Eq.(A.12), we find 
the following relation between C2 and Ci 

2'^^'^^ 3 
^2= 3(9^2 + i2/5«-m262)^i- ^^-^^^ 
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Using Eq.(A.17) reduces this to 



1 an 3 

which is an example of the general result 



and subsequently 



2n„2n 



Or, ~ K C 



(A.20) 



(A.21) 



(A.22) 



Thus all the coefficients are determined by the constant of integration ci. 
Solutions to the bulk equations must also satisfy the boundary conditions 



A' 



y=yi 



y=yi 



y=yi 



y=yi 



(A.23) 
(A.24) 



where i=l,2 refers to the boundary at y=0, and y=l respectively and Vi is 
specified below Eq.(A.l). We will first try to satisfy these boundary condi- 
tions by keeping only the first term in the series for $o and Ao and then use 
the relation Eq.(A.19) to determine if this is a reasonable approximation. 
Using 



(A.25) 
(A.26) 



in the $o boundary conditions we find equations for ci and e " in terms of 
vi and V2 



Cl 



rriiVi 



mi + ke 

7722^2 mi + ke 



777.2 — ke miVi 
Note that if 777i_2 ^ ke (the "stiff potential" limit) then 

Vl 



(A.27) 
(A.28) 

(A.29) 
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which was obtained in ^3] and used there to create a large hierarchy without 
the need for fine-tuning of parameters. Thus writing Eq.(A.29) as 



e-^''"^ (^)'^' (A.30) 

one see that for e.g. e = 1/30 one needs only assume V2/V1 = 0.3 to obtain 

g-feb. ^ 2 X 10~^^ (A.31) 

However, vi and V2 are not totally free parameters and using the boundary 
conditions for Ao we find relations for vi and f 2 



(A.32) 



2 6 nil + ek 
em I 

vl = 6 (A.33) 

67712 

When these relations are inserted into Eq.(A.28) we find that 

e-=f^V^^^V''. (A.34) 
ym2 — ke mi J 

This implies that e~" > 1 and therefore does not give a solution to the 
hierarchy problem. It should be noted that this result holds for any m,. 

One can also see that with e ^ 1 the n=2 term in the series is not small 
compared to the n=l term. From Eq.(A.19) the ratio C2/C1 becomes 

- ^ ^{^c^r (A.35) 
ci 12 



when e ^ 1. Using Eq.(A.27) for cf and Eq.(A.32) for vf we find 

C2 mi 
ci 2(mi + ke) 



(A.36) 



which is of 0{1), and in the stiff potential limit C2/C1 ~ 1/2. Therefore 
truncating the series to the first term is not a valid approximation, as the 
small parameter e in Eq.(A.35) cancels out in Eq.(A.36). It is the imposition 
of the A'^ boundary condition that makes Vi^2 ~ 1/e^/^. 
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We next examine the effect of retaining only the first two terms in the 
series expansions for and A^: 



Ao^ao + kbo + aie-'"^ + ase-^^^ 



(A.37) 
(A.38) 



As was previously noted, all coefficients can be found in terms of Ci. The $o 
boundary condition at y=0 becomes 



1 



l + 3(5i l + 35i 
where we have introduced the notation 

.2 \ 1/2 



vi = 



Ci = 



Vi = 



Si = 



12(1 + e); 

.2 \ 1/2 



Cl 



ere 



,12(1 + 



ke 



rrii 



(A.39) 

(A.40) 

(A.41) 
(A.42) 



Inserting Eqs.(A.37) and (A.38) into the y=0 boundary condition for A'^ gives 



SiCi (ci + Scfj = (ci + cf - Cl) . 



2(1 + e) 

Upon substituting for cf from Eq.(A.39), Eq.(A.43) becomes 
2(1 + S,)di - 3(1 - S,)c,v, + ^^^^^ - QS.vl = 0. 
This is easily solved for ci in terms of Vi. 



(A.44) 



Cl 



3vi 



4(l + (5i) 



l-5i± 



(3<5i + 1)' 



4(l + (5i)(l + 3(5i)- 
9{;?(1 + e) 



1/2N 



(A.45) 



In the stiff potential limit, 5i — > 0, this reduces to 



ci = ^|l± 



9vl{l + e) 



1/2N 



(A.46) 
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e 


0.03 


0.03 


0.01 


0.01 


Si 


0.01 


0.01 


0.0001 


0.0001 


52 


-0.01 


0.01 


-0.0001 


0.0001 


Vl 


±.6613 


±.6613 


±.6578 


±.6642 


V2 


±.6613 


± .6542 


±.6578 


±.6640 


Cl 


±.5152 


±.5152 


±.5184 


±.5219 




1 


0.9876 


1 


1.0001 




1 


0.6597 


1 


1.0101 



Table 1: Example of determination of the hierarchy parameter e~" for various 
choices of e, 6i, and 62- A vahd hierarchy is obtained when ~ 10~^^, which 
requires V1/V2 ~ 1/3 for e = .03 and Vi/v2 ~ 2/3 for e = .01 



Taking the positive root and putting this into Eq.(A.39) leads to an equation 
for Vl that can be solved to give 



-"1 



.667. 



(A.47) 



No real solution is found if one takes the negative root in Eq.(A.46). 
Prom the y=l boundary conditions we obtain the equations 



I-S2 
1-352 



1-3^2 



^2 = 



(A.48) 



3-;;2 



4(1 - ^2) 



1 + S2± 



(1 - 3S2f 



4(l-52)(l-352 



1 + e) 



1/2N 



(A.49) 

which are identical to Eqs.(A.39) and (A. 45) found at y=0 with ci replaced 
by e~"ci and 5i replaced by (— ^2)- In the stiff potential limit, S2 — > 0, this 
gives Vl — V2 and e"" = 1 which again would not give the desired solution 
to the hierarchy problem just as in the n=l case. 

We can also determine the situation for 6i small but non-zero. Table 1 
gives some sample values. Thus a hierarchy is not obtained if we truncate 
at n=2. The above results suggest that keeping a finite number of terms in 
Eqs.(A.lO) and (A. 11) will not lead to a valid approximation, and it may 
be that truncating at n=l does not approximate the rigorous solutions of 
Eqs.(A.7)-(A.9). 
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